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Fragmentation of a viscoelastic food by human mastication 
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Fragment-size distributions have been studied experimentally in masticated viscoelastic food 
(fish sausage). The mastication experiment in seven subjects was examined. We classified the 
obtained results into two groups, namely, a single lognormal distribution group and a lognormal 
distribution with exponential tail group. The facts suggest that the individual variability might 
affect the fragmentation pattern when the food sample has a much more complicated physical 
property. In particular, the latter result (lognormal distribution with exponential tail) indicates 
that the fragmentation pattern by human mastication for fish sausage is different from the 
fragmentation pattern for raw carrot shown in our previous study. 11 ' The excellent data fitting 
by the lognormal distribution with exponential tail implies that the fragmentation process has 
a size-segregation-structure between large and small parts. In order to explain this structure, 
we propose a mastication model for fish sausage based on stochastic processes. 
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1. Introduction 

Mastication is an in-mouth fragmentation process in 
which food is broken, ground or crushed by the teeth 
to prepare for swallowing and digestion. 1 ) In general, a 
major problem of the oral process observation analysis 
is lack of direct visualization of the process itself. This 
fact indicates that by investigating the food status before 
and after eating experimentally and numerically, we can 
understand some of the principal features of mastication 
processes and propose some phenomenological models. 2 ) 

As the first approximation, we propose that mastica- 
tion is a sequential fragmentation in the oral cavity be- 
tween the teeth and food. Fragmentation is a very com- 
plicated phenomenon, and it is difficult to understand its 
dynamics. One of methods used to understand the frag- 
mentation process is to investigate a fragment-size dis- 
tribution. The fragment-size distribution resulting from 
various types of fracture has attracted the interest of 
physicists for many years. We can give a few examples for 
the studies of the fragment-size distributions, the stud- 
ies of astelloids, 3 ) glass rods, 4 ) glass plates, 5 ) egg-shells 6 ' 
and so on. 

In dental science, there has been a considerable study 
of the fragment-size distribution. 7 ~ 10 ' We also experi- 
mentally studied the fragment-size distribution by mas- 
ticating raw carrots, which are regarded as a brittle ma- 
terial. 11 ' We reported that a lognormal distribution well 
fits the entire region for masticated fragments of raw car- 
rots for several chewing strokes. The above result indi- 
cates that the fragmentation of raw carrots by human 
mastication is characterized by the effect of random mul- 
tiplicative stochastic processes in statistics. 12 ' 

Solid foods which need mastication by the teeth to eat 
can be categorized into four groups by means of bite- force 



curve observed in the first chew of humans. 13 ' They are 
A) sponge-like foods (such as bread and sponge cake), B) 
gels (such as agar jelly, cooked rice, and fish gel), C) wet 
crisp foods (such as raw carrot, apple and cucumber), 
and D) dry crisp foods (such as cracker and cookie). As 
raw carrot belonging in the group C was studied in our 
previous paper, 11 ) we choose fish sausage from the gel 
group. Generally, gels are viscoelastic unlike crisp foods. 
The mechanical characteristics of viscoelastic gels were 
soft, easy to deform. Therefore, the viscoelastic food are 
even chewable for the elderly. 14 ' On the other hand, they 
are difficult to break off completely, 15, 16 - ) though it is ho- 
mogeneous and isotropic nature. We have a little infor- 
mation for the fragmentation process of viscoelastic food. 
The gel structure of fish sausage formed by fish protein 
molecular networks. 

Among the four groups discussed above, saliva absorp- 
tion effect during mastication may be greater in sponge- 
like food which has porous structure ready to hold liquid, 
and in dry-crisp food with high water absorption capac- 
ity. However since gel type food and wet-crisp foods are 
rich in water, no significant saliva absorption may oc- 
cur. In the model proposed by Hutchings and Lillford, 17 -' 
mastication is the process to break down the food struc- 
ture and to lubricate the bolus surface ready to swal- 
low. Evidently, the former is the main effect required for 
fish sausage and raw carrot. To discuss the differences in 
fragmentation by human mastication between those two 
groups of food is valuable. 

In this paper, we present an experimental study of 
the fragmentation of fish sausage by human mastication. 
So we examine data fittings by various distributions. Fi- 
nally, we discuss a plausible process of fragmentation of 
fish sausage by mastication in accordance with physical 
viewpoints. 
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2. Materials and Methods 

Test food and subjects. We used fish sausage 
(Osakanano-sausage, Nippon Suisan Kaisha, Lid., 
JAPAN) as test food in this study. The test food 
was cut into a cylinder (diameter and height; 24 mm 
and 15 mm, respectively) of about 7 g. Seven healthy 
subjects (4 males and 3 females, mean age 26.4 years) 
voluntarily participated in this study. They had natural 
dentition without severe malocclusion and periodontal 
disease. Written informed consent was obtained from 
each subject after full explanation of the experiment. 

Collect of food fragments. First, each subject mas- 
ticated the test food until swallowing as usual to count 
the number of chewing cycles until swallowing. Next, 
we calculated the individual number of chewing cycles 
until the halfway of the mastication. To collect the 
food bolus at the halfway of the mastication and just 
before swallowing, subjects were asked to masticate 
the test food until individually prescribed number of 
chewing cycles, and they spat the bolus into a beaker. 
The collection of food bolus was performed two times 
per subject and each condition. In order to expectorate 
entirely the subject rinsed their mouth with water. The 
food fragments and water were carefully stirred in a 
beaker with a grass rod and pass through a sieve with a 
mesh size of 0.5 mm. After fine fragments were washed 
through the sieve with running water, the fragments on 
the sieve were spread evenly on the transparent acrylic 
board (300 mm x 300 mm). Then, we made a copy of 
this board using copy machine (DocuCentra-II, FUJI 
XEROX, JAPAN) without cover and stored the copy on 
a personal computer at a suitable resolution (about 0.13 
mm/pixel). 

Statistics. The additional data from second trials 
was then assimilated and sorted in terms of size from 
the largest one, i.e., cumulative number. Then, we ob- 
tained the cumulative size distribution for each number 
of chewing strokes and applied a suitable distribution to 
fit a curve to the data in each case. For fitting the curve, 
we performed a nonlinear least-square method using R 
(version 2.4.1 for windows). 1S ) 

3. Results and Discussion 

Figures 1-7 show the cumulative number of food- 
fragments. It should be noted that different figures cor- 
respond to different subjects, respectively. These figures 
are classified roughly into two groups based on fitting 
curves, namely, single lognormal distribution (Figs. 1, 2, 
6(a) and 7(a)) and lognormal distribution with an expo- 
nential tail (Figs. 3, 4, 5, 6(b) and 7(b)). For example, as 
Fig. 1(a) shows, the fragment-size distribution is nicely 
approximated by a single lognormal distribution given as 
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Fig. 1. Log-log plots for the cumulative numberof masticated 
food fragments of fish sausage after (a) 32 and (b) 16 chews. 
The solid line indecates a lognormal distribution for (a) = 
1161,5 = 0.0375,ct = 0.859 and N/N T = 1.08, and (b) N T = 
688, 5 = 0.0559, o = 0.965 and N/N T = 1.04. 
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Fig. 2. Same as Fig. 1 except for the number of chewing strokes, 
i.e., (a) 24 and (b) 12. The fitting parameters are (a) JVt = 
620,5 = 0.0535,cr = 1.10 and N/N T = 1.02, and (b) N T = 
319, s = 0.0837, o = 1.13 and N/N T = 1.02. 
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Fig. 3. Same as Fig. 1 except for the numner of chewing strokes, 
i.e., (a) 50 and (b) 25. The fitting parameters are (a) JVt = 
1162,5 = 0.0341,o = 0.754 and N/N T = 1.11, and (b) N T = 
784,5 = 0.0369,o- = 1.11 and N/N T = 1.15. The insets show 
semi-log plots for the tail part. The fitting parameters are (a) 
A = 0.732, B = 0.0497 and (b) A = 0.474, B = 0.118. 



where a and s are the fitting parameters, respectively. 
The cumulative form for the lognormal distribution is 



N(s) = -^(1 



e rf( 1 ° g( / f /s " ) )), 



(2) 



V2a 

where Nt is the total number of fragments and 
erf(x) is the error function defined as erf(x) = 
(2/-^) J* exp(—y 2 )dy. A similar result was reported in 
our previous study using raw carrot. 11/1 

On the other hand, as shown in Figs. 3, the majority 
of distributions belong to the small region approximated 
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(a) (b) 

Fig. 4. Same as Fig. 1 except for the number of chewing strokes, 
i.e., (a) 70 and (b) 35. The fitting parameters are (a) = 
1162, 5 = 0.0260,ct = 0.750 and N/N T = 1.12, and (b) N T = 
873,5 = 0.0402,ct = 0.964 and N/N T = 1.16. The fitting pa- 
rameters for the insets are (a) A = 0.568, B = 0.0425 and (b) 
A = 0.531, B = 0.103. 
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Fig. 5. Same as Fig. 1 except for the number of chewing strokes, 
i.e., (a) 58 and (b) 29. The fitting parameters are (a) Ni = 
1167, s = 0.0355,o = 0.729 and N/N T = 1.08, and (b) N T = 
611,5 = 0.0471,o = 1.25 and N/N T = 1.15. The fitting pa- 
rameters for the insets are (a) A = 0.648, B = 0.0562 and (b) 
A = 0.643, B = 0.139. 
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Fig. 6. Same as Fig. 1 except for the number of chewing strokes, 
i.e., (a) 30 and (b) 15. The fitting parameters are (a) = 
576,5 = 0.0595,o- = 1.03 and N/N T = 1.03, and (b) N T = 
371,5 = 0.0384,o- = 1.55 and N/N T = 1.30. The fitting param- 
eters for the inset arc (b) A = 0.468, B = 0.248. 

by the lognormal distribution. However, the tail part of 
the distributions deviates upwards from the lognormal 
distribution. Hence, we propose that the size segregation 
into small and large food fragment groups occured due 
to physical or other properties of fish sausage. The insets 
of Figs. 3 show semi- log plots for the tail part (same as 
Figs. 4, 5, 6(b) and 7(b)). Since the curves on the insets 
are plotted linearly, the large food fragments were fitted 
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Fig. 7. Same as Fig. 1 except for the number of chewing strokes, 
i.e., (a) 20 and (b) 10. The fitting parameters are (a) JVt = 
722,5 = 0.0468,o- = 1.06 and N/N T = 1.15, and (b) N T = 
326, 5 = 0.0444,0- = 1.34 and N/N T = 1.34. The fitting param- 
eters for the inset arc (b) A = 0.439, B = 0.236. 

to not a lognormal but an exponential distribution, 

N(s) = Ae-^, (3) 

where A and B are the fitting parameters, respectively. 
The physical origin of an exponential distribution of frag- 
mentation is very simple. If we assume that a fragment- 
size at each stage of sequential fragmentation is com- 
pletely random, then we obtain the exponential distri- 
bution as the fragment-size distribution. 19 ) In the case 
of our mastication experiments, we think that a similar 
phenomenon like the above assumption happens. 

Using the above results, we propose a mechanism for 
the mastication process of fish sausages (see Fig. 8). In 
the first stage of mastication, a unit of fish sausage is 
broken into a small amount of fragments. Fish sausage 
is then divided into several more fragments without gen- 
erating small fragments, while raw carrot is divided into 
many fragments of varied sizes. 11 ) The difference of the 
fragmentation pattern between fish sausage and raw car- 
rot is derived from their physical properties that is fish 
sausage behaves more plastic than raw carrot in the 
mouth If the generated fragments of fish sausage are 
larger than about 0.1 cm 2 , the fragments are divided into 
several fragments again (exponential in Fig. 8). If the gen- 
erated fragments are smaller, the fragments are ground 
by the back teeth, and then many, heterogeneous and 
smaller fragments are generated. After several cycles, the 
fragments caused by the successive mastication process 
for fish sausage will satisfy the two threshold hypothesis 
suggested by J. B. Hutchings and P. J. Lillford. 17 ) Here, 
the chewed food is assembled into a bolus by a compli- 
cated movement of the palate and the tongue, then just 
before swallowing. 20 ) According to this proposition, we 
will conjecture that a fragment-size distribution obtained 
by grinding fish sausage leads to a lognormal distribu- 
tion (lognormal in Fig. 8). In fact, Epstein theoretically 
showed that the fragment-size distribution by grinding 
solid leads to lognormal. 21 ' The proof of this conjecture 
is one of our future problems. 

As mentioned above, we conclude that there are two 
different kinds of fragment-size distributions, i.e., single 
lognormal distribution and lognormal distribution with 
exponential tail, because the tail part of distribution 
drifted either downwards or upwards from a lognormal 
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Fig. 8. A flow chart of the mastication process for fish sausages. 



distribution fitted to a small fragment group. This fact 
indicates that there are differences, such as saliva produc- 
tion and so on, in the mastication process for fish sausage 
among individuals, unlike for raw carrot. 11 ) In fact, indi- 
vidual subjects have shown different behaviour in man- 
aging food to terminal swallow. 20 > The second possibility 
is that the single lognormal distribution could be close 
to the lognormal distribution with exponential tail in re- 
peated experiments because we must discuss about the 
probability density distribution. This is also our second 
future problem. 

4. Conclusion 

In summary, we have studied the fragment-size dis- 
tribution of masticated viscoelastic food (fish sausage). 
We classified the obtained results from 7 humans into 
two groups, namely, the single lognormal group and the 
lognormal distribution with exponential tail. The former 
group shows a similar result as the one in our previous 
studies using raw carrot. The latter group is the particu- 
lar case of masticated fish sausage. The fragment-size dis- 
tribution for the latter group shows a double-size-group 
tendency, i.e., the majority of the distribution belongs 
to the lognormal distribution, while the tail behaves as 
the exponential distribution. In order to explain this ten- 
dency, we suggest a mastication model for fish sausage 
shown in Fig. 8. However, the mastication model remains 



unfinished, and then we have a lot of further problems, 
such as the size distribution by grinding fragmentation 
among others. 
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